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ABSTRACT 

We analyze deep g' and r' band data of 97 galaxy clusters imaged with McgaCam on 
the Canada- France-Hawaii telescope. We compute the number of luminous (giant) and 
faint (dwarf) galaxies using criteria based on the definitions of de Lucia et al. (2007). 
Due to excellent image quality and uniformity of the data and analysis, we probe the 
giant-to-dwarf ratio (GDR) out to z ~ 0.55. With X-ray temperature (T x ) information 
for the majority of our clusters, we constrain, for the first time, the T^-corrected giant 
and dwarf evolution separately. Our measurements support an evolving GDR over the 
rcdshift range 0.05 ^ z ^ 0.55. We show that modifying the (g' — r'), m r * and K- 
correction used to define dwarf and giant selection do not alter the conclusion regarding 
the presence of evolution. We parameterize the GDR evolution using a linear function 
of redshift (GDR = az + (3) with a best fit slope of a — 0.88 ± 0.15 and normalization 
/3 = 0.44 ± 0.03. Contrary to claims of a large intrinsic scatter (<7mt), we find that the 
GDR data can be fully accounted for using observational errors alone. Consistently, 
we find no evidence for a correlation between GDR and cluster mass (via T x or weak 
lensing). Lastly, the data suggest that the evolution of the GDR at z < 0.2 is driven 
primarily by dry merging of the massive giant galaxies, which when considered with 
previous results at higher redshift, suggests a change in the dominant mechanism that 
mediates the GDR. 

Key words: galaxies: clusters: general - galaxies: evolution - galaxies: elliptical and 
lenticular, cD 



1 INTRODUCTION 

Non-starforming galaxies in clusters exhibit a tight sequence 
in colour-magnitude space known as the red sequence (eg. 
Visvanathan 1978; Bower et al. 1992), which is in place 
as early as z ~ 1 (eg. Bell et al. 2004, Mei et al. 2006). 
Less understood is how the number of red sequence galaxies 
and their distribution with stellar mass/absolute magnitude 
evolves with cosmic time. Understanding the assembly his- 
tory of the red sequence is critical as these galaxies represent 
the bulk of the stellar mass in clusters at the present epoch. 
A useful metric for probing this assembly history is the ratio 
of the number of luminous red sequence galaxies (giants) to 
the number of faint red sequence galaxies (dwarfs). This is 
essentially a non-parametric representation of the luminos- 
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ity function using only 2 bins and is commonly referred to 
as the Giant-to- Dwarf ratio (GDR). 

There is an ongoing debate regarding the presence or 
absence of evolution in the GDR and/or its reciprocal the 
Dwarf-to-Giant ratio (DGR). De Lucia et al. (2007) for in- 
stance, finds significant evolution in the GDR, amounting 
to a change from GDR~ 0.45 to GDR~ 0.95 over the red- 
shift range 0.4 < z < 0.8. The analysis of low redshift 
(0.08 < 2 < 0.19) cluster data by Capozzi et al. (2010) sup- 
ports the extrapolation of this trend to lower redshift. Stott 
et al. (2007) and Gilbank et al. (2008) also find strong evo- 
lution in the DGR using somewhat brighter absolute mag- 
nitude limits for the definition of dwarfs. In contrast, Craw- 
ford et al. (2009) look at 59 clusters at redshift < z < 0.5 
and, though they do not measure the GDR directly, con- 
clude that there is little evolution in the faint-end slope of 
the luminosity function. Andreon (2008), in an analysis of 
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28 clusters in the redshift range < z < 1.3, asserts that 
the DGR is consistent with no evolution. More specifically, 
he states that the strong trend observed by de Lucia et al. 
(2007) can be ruled out. Andrcon (2008) discusses some of 
the differences between his analysis and those of de Lucia et 
al. (2007) and Gilbank et al. (2008), arguing that their selec- 
tion of filters that do not bracket the 4000 angstrom break 
over the entire range of redshifts introduces systematic er- 
rors that mimic evolution. Furthermore, he argues that in 
the work of Stott et al. (2007) and de Lucia et al. (2007) 
the use of control fields taken in separate filter sets than the 
cluster fields causes significant systematics related to the 
interloper-removal process. 

The selection of target clusters may also be important. 
The Andreon (2008) clusters are all spectroscopically con- 
firmed and identified in X-rays, while those used in de Lucia 
et al. (2007) are not and may contain line-of-sight superpo- 
sitions at the low-mass end. 

Perhaps the most significant difference between An- 
dreon (2008) and others is the use of a Bayesian approach to 
fit the luminosity function of red sequence galaxies, which is 
then converted to a DGR to facilitate comparison with the 
literature. It is stated that this is a superior method because 
it uses all of the available data, particularly at the faint-end 
of the luminosity function, beyond the limits of the de Lucia 
et al. (2007) definition of a dwarf galaxy (the most common 
definition). We suspect however, that there are distinct dis- 
advantages to the Andreon (2008) type approach because 
in forcing the parametric form of the luminosity function 
one is introducing an undesirable model dependence and it 
may therefore be incorrect to compare with non-parametric 
results. Moreover, a procedure that takes all available data 
and allows the fitted domain to vary from cluster to clus- 
ter potentially creates a systematic bias. This is because the 
clusters that are sampled to the faintest absolute magni- 
tudes will have more stringent constraints on their faint end 
slope than those that are truncated at a brighter limit. 

With so many subtle differences between the various 
ways that the GDR has been measured in the literature 
it may be no surprise that there appears to be some dis- 
agreement. Our analysis, however, is carried out on ob- 
servations taken in a single filter pair, from a single tele- 
scope/instrument (including control fields), and with a sam- 
ple of exclusively X-ray-confirmed clusters. These attributes 
allow us to avoid many of the criticisms mentioned above. 
Moreover, we calculate the effect of many potential sources 
of systematics (7f -correction, projected distance from the 
brightest cluster galaxy (BCG), colour selection method, 
colour error bias and the presence of large scale structure) 
to estimate their role in the apparent disagreements in the 
literature. When coupled with the sheer number of clusters 
in our sample, this survey provides the most accurate and 
robust measure of the GDR over this redshift range to date. 

If indeed there is an evolving GDR in galaxy clusters, it 
becomes interesting to investigate the mechanism(s) respon- 
sible for driving that evolution. For instance, gas-poor (dry) 
merging of galaxies at the bright end of the red sequence 
could lead to a depletion of the giant population. In compe- 
tition with this process, star-forming galaxies move onto the 
red sequence through the quenching of ongoing star forma- 
tion and the resulting colour transformation. However, the 
observed downsizing behavior of stellar mass assembly (Bell 



et al. 2004, Juneau et al. 2005, Bundy et al. 2005, Faber 
et al. 2007, Scarlata et al. 2007) suggests that this effect 
may become less important with decreasing redshift as the 
specific star formation rate at a fixed galaxy mass declines. 
While these mechanisms are complementary in that they act 
to reduce the number of giants with respect to dwarfs, their 
relative importance as a function of look-back time is poorly 
constrained. The detailed behaviour of these mechanisms is 
likely an important benchmark for galaxy formation models. 
For example, Gilbank & Balogh (2008) show that the Gal- 
form model of Bower et al. (2006) is unable to reproduce 
both the observed DGR evolution and the DGR difference 
between cluster and field populations. 

As mentioned above, all of our clusters are confirmed 
in X-rays and the majority of them have X-ray tempera- 
ture measurements, which can be used as independent mass 
proxies. This information is necessary for disentangling any 
cluster mass dependence and allows us, for the first time, to 
properly probe the mechanism(s) that drive GDR evolution 
at these redshifts. 

The layout of the paper is as follows. In f[5]we discuss 
the data, along with some tests for completeness and data 
quality to ensure that our measurements are not significantly 
biased by systematics. In §3] we describe the data reduction 
and analysis procedure used to measure the red sequence 
and extract from it the GDR. Following this, i|4] discusses 
the results of our GDR analysis along with tests of how 
varying the selection method can affect the GDR measure- 
ment and a detailed comparison with previous results from 
the literature. In iJS]we discuss the potential dependence of 
the GDR on cluster mass as implied by the X-ray temper- 
ature. An investigation of the mechanisms responsible for 
GDR evolution is in SjfJ] Finally, iJ7]contains our conclusions. 
Throughout this paper we assume a cosmology parameter- 
ized by f2 m = 0.3, A = 0.7 and H = 70km/s/Mpc. 



2 DATA 

The data for this paper are obtained using the MegaCanQ 
wide-field imager on the Canada-France-Hawaii Telescope 
(CFHT) with the g' and r' optical filter set. The cluster 
sample is comprised of 101 clusters with redshift 0.05 < 
z < 0.55 and is assembled from 3 subsamples: 58 clusters 
from the Multi-Epoch Nearby Cluster Survey (MENeaCS), 
32 clusters from the Canadian Cluster Comparison Project 
(CCCP) and 11 additional clusters from the CFHT archive. 
We remove 4 clusters from the analysis because of their low 
richness (see £|4.1fl leaving a total of 97. Basic information 
for our targets is listed in Table Q] Further details on the 
CCCP survey will be presented in a future paper (Hoekstra 
in prep.), while a more detailed description of the MENeaCS 
survey can be found in Sand et al. (2012). The largest seeing 
disc of our images has a FWHM of 1.0" with the majority be- 
ing significantly smaller. For background determination, we 

1 Based on observations obtained with MegaPrimc/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institut National 
des Science de l'Univers of the Centre National de la Recherche 
Scicntifique (CNRS) of France, and the University of Hawaii. 
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Figure 1. Half-light radius vs apparent magnitude in the r' band 
for the clusters Abcllll9 at z = 0.05(a) and RXJ1524 at z = 
0.516 (b & c). The black polygons illustrate our star removal 
procedure. The red triangles and blue squares represent all of 
the giant and dwarf galaxies respectively (see i|4] for details of 
the selection procedure). For clarity we show only 5% of regular 
objects but 100% of the selected giants and dwarfs. The m r i giant 
and dwarf selection limits are shown by the red and blue dashed 
lines. The solid vertical magenta line shows the equivalent size 
of a z = 0.05 galaxy measuring 10" after taking into account 
the PSF and the change in angular diameter distance. The (c) 
panel shows only the giants, dwarfs and fainter galaxies (black 
pentagons) that also matched the red sequence colour selection 
(see SQJ. The large size of these fainter galaxies with respect to 
the stellar distribution indicates a clear separation of stars and 
red sequence galaxies based on our selection criteria. 



make use of the CFHT Legacy Survey (CFHTLS) Deep im- 
ages. These data are processed and hosted by S. GwyrQ (see 
Gwyn 2008 for details) and consist of four fields observed 
away from any galaxy clusters. The CFHTLS Deep fields are 
observed using the same instrument and filters as our tar- 
geted observations but to greater depth. As such, they pro- 
vide excellent control fields for statistical background sub- 
traction (discussed further in Q. 



2.1 Source Extraction 

The data are bias-subtracted and flat fielded using the Elixir 
software package. We combine the images and perform back- 
ground subtraction using SWarp (Bertin et al. 2002). After 
data reduction and image stacking we mask obvious image 
defects and regions around bright stars that could severely 
bias the photometry. We then run SExtractor (Bertin & 
Arnouts 1996) in dual image mode on each set of g' and 



2 http:/ /wwwl.cadc-ccda.hia-iha.nrc- 

cnrc.gc.ca/community/CFHTLS-SG/ docs / cfhtlsDccp.html 




1TL , 



Figure 2. Histograms showing the surface-density of galaxies per 
sq. arcmin for 4 of the highest redshift clusters in our sample. The 
thick lines in each panel show the galaxies detected within a wide 
(red) and narrow (blue) colour slice (see text). The correspond- 
ing thin black lines show the surface-density of galaxies detected 
in the CFHTLS Deep images using the same selection criteria. 
The blue vertical dashed line shows the location of the faint-end 
magnitude limit used to select dwarf galaxies (see j|4]for details). 
The data are complete up to this limit for the full redshift range 
investigated. 



r' images using the masks merged with exposure maps as 
weight images. It is critical that we detect and photometer 
all faint dwarf galaxies, even when they are embedded in the 
halos of the much larger giant galaxies and so we adjust some 
of the SExtractor parameters from their default settings ac- 
cordingly. We briefly discuss here some of the SExtractor 
parameters that can affect the number of sources detected 
in crowded regions and the values we use for our extractions. 
The DETECT.THRESH and ANALYSIS.THRESH values 
are set to 1.5 relative to the background variance and the 
DETECT_MINAREA is set to 3 pixels. To enable the detec- 
tion of faint galaxies near to bright ones in these crowded 
fields we set the DEBLEND_MINCONT to 10~ 7 and the 
DEBLEND JNfTHRESH to 64. There are also several param- 
eters that control the way the background is measured and 
subtracted from the detected sources and, therefore, these 
can affect the apparent magnitudes of galaxies especially at 
the faint end. For these we use a BACK_SIZE of 128 pixels 
and a BACK_FILTERSIZE of 3 pixels. Lastly, SExtractor 
has the option to "CLEAN" the source catalog of artifacts 
and false detections. We find that setting CLEAN =Y with 
a CLEAN_PARAM= 2.0 not only removes many of the ar- 
tifacts and false detections but also increases the number of 
real objects detected compared to the case of CLEAN=N. 

Throughout this paper we use SExtractor MAG_AUTO 
for galaxy magnitudes and MAG_APER for galaxy colours. 
The aperture magnitudes are extracted using a 3" diameter 
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Table 1. Table of clusters used in our sample. RA (a) and DEC (8 
* are excluded from the final analysis due to low richness (sec JJTJ 



Name z a 8 

hh:mm:ss dd:mm:ss 



AbeM19 





.044 


00 


56 


16 


,10 


-01 


15 


19 





MKW3S 





,045 


15 


21 


51 


,81 


07 


12 


31 


9 


Abell780 





,054 


09 


18 


05 


65 


-12 


,05 


13 


.5 


Abell754 





,054 


09 


08 


32 


37 


-09 


37 


47 


2 


Abell85 





,055 


00 


41 


50 


,45 


-09 


:18: 


11 


,1 


Abell2319 





,056 


19 


21 


10 


00 


13 


,56 


11 


.5 


Abelll33 





,057 


01 


02 


41 


71 


-21 


52 


55 


2 


Abelll991 





,059 


14 


51 


31 


,48 


18 


38 


33 


4 


Abelll781* 





,062 


13 


:44 


52 


,54 


29 


46 


15 


,6 


Abelll795 





,062 


13 


18 


52 


49 


26 


35 


31 


,8 


Abell553 





,066 


06 


12 


11 


09 


18 


35 


11 


6 


Abell644 





,070 


08 


:17 


25 


61 


-07 


30 


45 





Abell399 





,072 


02 


:57 


53 


,09 


13 


01 


51 


3 


Abell2065 





,073 


15 


22 


29 


,17 


27 


12 


27 


,8 


Abell401 





,074 


02 


,58 


57 


,78 


13 


31 


58 


3 


ZwC11215 





,075 


12 


:17 


41 


13 


03 


39 


21 


2 


Abell2670 





,076 


23 


51 


13 


,67 


-10 


25 


08 


2 


Abell2029 





,077 


15 


10 


56 


09 


05 


11 


11 


4 


Abell2495 





,078 


22 


50 


19 


,71 


10 


51 


12 


,8 


RXSJ2344-04 





,079 


23 


:44 


18 


20 


-01 


22 


18 


,8 


ZwC10628 





,081 


06 


31 


22 


69 


25 


01 


06 


,8 


Abell2033 





,082 


15 


11 


26 


51 


06 


20 


56 


,8 


Abelll650 





,084 


12 


,58 


41 


49 


-01 


15 


11 


,0 


AbehT651 





,085 


12 


,59 


22 


49 


-01 


11 


15 


,7 


Abell2420 





,085 


22 


10 


18 


,76 


-12 


10 


13 


9 


Abell2597 





,085 


23 


:25 


19 


72 


-12 


:07 


26 


,7 


Abell763* 





,085 


09 


12 


35 


,18 


16 


00 


01 


,0 


Abell478 





,088 


01 


;13: 


25 


27 


10 


27 


55 


,1 


Abell2440 





,091 


22 


23 


56 


92 


-01 


34 


59 


.5 


Abell2142 





,091 


15 


:58 


19 


99 


27 


;14: 


00 


.5 


Abelll927 





,095 


11 


31 


06 


,79 


25 


38 


01 


6 


Abell21 





,095 


00 


20 


36 


98 


28 


39 


33 





Abell2426 





,098 


22 


11 


31 


,58 


-10: 


22 


26 


,1 


Abell2055 





,102 


15 


:18: 


45 


72 


06 


13: 


56 


4 


Abelll285 





106 


11 


30 


23 


,82 


-14 


31 


52 


3 


Abell7 





106 


00 


11 


15 


25 


32 


21 


56 


6 


Abell2064* 





,108 


15 


20 


52 


,24 


48 


39 


38 


,7 


Abell2443 





,108 


22 


26 


07 


92 


17 


21 


23 


,8 


RXCJ0352+19 





,109 


03 


:52 


:58 


.99 


19 


10 


:59 


,8 


Abell2703 





,114 


00 


05 


23 


95 


16 


:13: 


09 


3 


Abell2069 





,116 


15 


21 


08 


42 


29 


52 


55 


,6 


AbeU1361 





117 


11 


43 


39 


60 


46 


21 


20 


,7 


Abell2050 





118 


15 


16 


17 


92 


00 


,05 


20 


9 


RXCJ0736+39 





118 


07 


:36 


:38, 


,08 


39 


:24 


:52 


,8 


Abelll348 





119 


11 


41 


21 


19 


-12 


16 


38 


.5 


Abell961 





,124 


10 


16 


22 


,86 


33 


38 


17 


,7 


Abelll033 





126 


10 


31 


11 


32 


35 


02 


29 


2 


Abell2627 





126 


23 


36 


12 


,08 


23 


,55 


29 


.5 


Abell655 





,127 


08 


:25 


29 


,05 


47 


08 


00 


9 



aperture, which is several times larger than the worst case 
seeing disc of the sample (see Sj2j. All catalogs are corrected 
for foreground galactic extinction using E(B-V) values from 
Schlegel et al. (1998) and the extinction curve of Cardelli et 
al. (1989). 



give the coordinates of the BCG. Cluster identifiers marked with a 



Name z a 8 

hh:mm:ss dd:mm:ss 



Abell646 





129 


08: 


22 


09 


,53 


47: 


05 


:53 


3 


Abellll32 





136 


10 


58 


23 


65 


56 


17 


12 





Abell795 





136 


09 


21 


,05 


29 


14: 


10 


:21 


,8 


Abelll068 





,138 


10 


40 


11 


48 


39 


57 


11 


,5 


Abelll413 





,143 


11 


55 


:18, 


00 


23 


21 


:18 


,1 


ZwC11023 





,143 


10 


25 


,57 


98 


12 


11 


,08 


7 


Abell990 





144 


10 


23 


39 


91 


19 


08 


:38 


,8 


Abell2409 





,148 


22 


00 


53 


49 


20 


58 


: 12 


,1 


RXCJ0132-08* 





,149 


01 


32 


41 


11 


-08: 


01 


01 


6 


Abell2204 





150 


16 


32 


46 


96 


05 


31 


33 


.0 


Abell545 





154 


05 


32 


25 


18 


-11 


32 


39 


3 


Abell2104 





,160 


15 


40 


:07 


92 


-03 


18 


16 


,1 


Abell2259 





,160 


17 


20 


09 


65 


27 


40 


:08 


,5 


Abelll234 





,166 


11 


22 


29 


,94 


21 


21 


22 


,2 


Abelll914 





170 


14: 


26 


03 


,88 


37 


19 


53 


,8 


Abell586 





,170 


07 


32 


20 


30 


31 


38 


:01 


3 


Abelll246 





190 


11 


23 


:58 


82 


21 


28 


50 


,0 


MS0440 





190 


01 


13 


09 


90 


02 


10 


19 


.5 


Abellll5S 





,200 


00 


56 


00 


25 


26 


20 


32 


,8 


Abellll5N 





,200 


00 


55 


50 


61 


26 


21 


37 


.7 


Abell2163 





,200 


16 


15 


:48 


98 


-06: 


08 


41 





Abell2261 





,200 


17 


22 


:27 


22 


32 


07 


:57 


,9 


Abell520 





,200 


01 


51 


:14, 


05 


02 


57 


10 


,8 


Abell223N 





,210 


01 


38 


02 


28 


-12 


,15 


19 


7 


Abell223S 





,210 


01 


37 


:55 


98 


-12: 


49 


09 


9 


Abelll942 





,220 


14 


38 


21 


87 


03 


40 


13 


,5 


Abell2111 





,230 


15 


39 


40 


50 


31 


25 


,27 


7 


Abell2125 





,246 


15 


11 


:14, 


75 


66 


16 


03 


,8 


Abelll835 





,250 


14 


01 


02 


09 


02 


52 


■12 


,9 


Abell521 





,250 


04 


51 


06 


87 


-10: 


13 


2-1 


4 


CL1938 





,260 


19 


38 


18 


10 


51 


09 


40 


,4 


Abelll758W 





,280 


13: 


32 


:38 


11 


50 


:33 


:36 


,2 


Abell697 





,280 


08 


42 


:57 


56 


36 


21 


:59 


6 


Abell611 





,290 


08 


00 


56 


82 


36 


03 


:24 





Abell959 





,290 


10 


17 


36 


00 


59 


34 


:02 


.0 


Abell2537 





,300 


23 


08 


22 


21 


-02: 


11 


31 


5 


MS1008 





,300 


10 


10 


32 


32 


-12: 


39 


52 


6 


Abell851 





,410 


09 


42 


:57 


45 


46 


58 


:50 


1 


RXJ0856 





,411 


08 


56 


12 


69 


37 


56 


15 


9 


RXJ2228 





,412 


22 


28 


33 


70 


20 


37 


16 


6 


RXJ1003 





,416 


10 


03 


01 


62 


32 


53 


41 


4 


MS1621 





,426 


16 


23 


35 


14 


26 


34 


:28 


3 


MACS1206 





,440 


12 


06 


12 


,14 


-08 


18 


03 


1 


CL0910 





,440 


09 


13 


15 


50 


10 


56 


,28 


7 


RXJ1347 





,450 


13 


47 


30 


64 


-11: 


45 


08 


9 


RXJ1701 





,453 


17 


01 


23 


51 


64 


11 


12 





3C295 





,460 


14: 


11 


20 


55 


52 


12 


10 


1 


RXJ1524 





,516 


15 


21 


38 


37 


09 


57 


43 


6 


MS0451 





,540 


04 


51 


10 


83 


-03: 


00 


:51 


2 



2.2 Stars & Extended Galaxies 

The object catalogs discussed in §2. II include stars. To ob- 
tain a galaxy catalog we remove the objects that lie along 
the locus defined by point sources in the half-light radius 
(?"o.5) vs apparent r' magnitude (m r /) plane. Figure [T] illus- 
trates this process for two clusters in our sample; one at the 
lowest redshift and one at the highest. The thick, solid, black 
box outlines the objects that are identifiable as unsaturated 
stars while the thin dotted polygon similarly identifies the 
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Figure 3. The average observed (g' — r') colour of red sequence 
galaxies with —20.82 < M r i < —19.82 determined using a boot- 
strap method. This absolute magnitude bin corresponds to half a 
mag on either side of the division between dwarfs and giants. The 
data are plotted as a function of redshift. The strong correlation 
observed is the effect of the g' filter moving across the 4000 A 
break with increasing redshift. 
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Figure 4. The average K-corrected (g' — r') colour of red se- 
quence galaxies with —20.82 < M r i < —19.82 determined using a 
bootstrap method. The data are plotted as a function of redshift. 
It is clear that after i\-correction the mean (g' — r') colour at the 
division between dwarf and giant galaxies is roughly constant. 



saturated ones. At magnitudes fainter than the limit of the 
solid box it becomes impossible to confidently distinguish 
point sources from extended ones by this method. For this 
reason, we retain all sources fainter than this limit in our 
galaxy catalogs so as not to remove the faintest galaxies. 
We verify, however, that any remaining stellar contamina- 
tion does not bias our results by examining whether the 
GDR (as measured in Q changes when excluding clusters 
at galactic latitude b < 35 degrees, which it does not. 

The blue squares and red triangles in Figure [T] show 
the red sequence dwarf and giant galaxies respectively, with 
the blue and red dashed lines showing the correspond- 
ing magnitude limits for their selection. These limits are 
adopted from the de Lucia et al. (2007) definitions (dwarfs: 
-20.0 < M v < -18.2 & giants: M v < -20.0) which are 
converted to the corresponding m r i values at the appropri- 
ate redshift (see lj4]for details of this conversion). We note 
the presence of a small number of detections with ro.s > 1.0" 
and 77V > 24.0 for some images (eg. RXJ1524). We have in- 
vestigated these by eye and conclude that they are all false 
detections and, because they do not enter the dwarf and 
giant selection limits, they do not affect the GDR. 

To verify that the number of dwarfs or giants is not 
affected by the detection efficiency of galaxies with rela- 
tively low central surface brightness, we examine the most 
extended galaxies near the detection limit in each cluster. 
From visual inspection of ro.s vs m r i (eg. Figure[T}, we find 
that at z ~ 0.05 there are very few galaxies detected that are 
larger than 10". At z ~ 0.05, after subtracting in quadra- 
ture the size of the PSF, we estimate the physical size of 
the limiting half-light radius to be ro.s ~ 8.5 kpc. This limit 
is then converted back to an angular size at the appropri- 
ate redshift, corrected for PSF effects and then overplotted 
as the solid magenta line in Figure [1] The analogs of these 
extended galaxies are well above our detection limit in the 
high z end of our sample. We conclude that the evolution 
of the number of dwarf and giant galaxies is not affected by 
the detectability of dwarfs with extended surface brightness 
profiles. Furthermore the fraction of such galaxies is small 
in both the dwarf and giant populations. 

2.3 Completeness 

To investigate the completeness of our galaxy catalogs we 
compare them with the catalogs obtained from the CFHTLS 
Deep data (see Ej2j) . We generate Galactic extinction- 
corrected, star-cleaned catalogs in both g' and r' for each 
of the 4 Deep fields using the same procedure as described 
in sections 12.11 and 12.21 We compare the number of galax- 
ies as a function of magnitude m r i in our sample with 
those in the CFHTLS Deep (control) using both a wide 
colour window [—1.0 < (g' — r') < 4.0] and a narrow one 
[1.0 < (<?' — r') < 1.2]. Figure [2] shows the surface den- 
sity of galaxies selected with these criteria for 4 clusters at 
the upper redshift limit of our sample. The thick red and 
blue histograms show the distributions of sources outside 
a projected radius of 750 kpc in each of our cluster fields 
for the wide colour window (red) and the narrow window 
(blue) while the underlying thin black histograms show the 
corresponding distributions of sources in the CFHTLS Deep 
survey area. The vertical blue dashed lines show the loca- 
tion of the faint end limit for dwarf galaxies (see The 
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slopes of the counts distributions are in agreement between 
our cluster data and the CFHTLS Deep fields. Some clus- 
ter fields show a small excess of objects with respect to the 
control sample, but this is expected due to enhancement 
from cluster galaxies that lie further than 750 kpc from the 
BCG. Some fields also exhibit a small jump at m T i ~ 22 
(eg. RXJ1524, MS0451) which corresponds to the apparent 
magnitude where we can no longer reliably distinguish stars 
from galaxies. The clusters that exhibit the largest jumps 
are those in the galactic plane b < 35°. As mentioned in 
i|2.2l we verified that these clusters do not bias our results 
by repeating our analysis with b < 35° removed and ob- 
taining consistent measurements of the GDR evolution. The 
CFHTLS Deep data is several magnitudes deeper in the r' 
band than our cluster sample, thus the surface- density be- 
gins to fall off in the clusters fields at brighter m r > than in 
the control fields. In all cases however, the sharp drop in 
counts observed in our sample occurs at magnitudes fainter 
than the faint dwarf limit (vertical dashed line). We con- 
clude that the GDR measurements presented in iJH are not 
affected by incompleteness in the dwarf population. 



3 ANALYSIS 

In this section we describe the various criteria we use to 
select galaxies belonging to the dwarf and giant populations 
within clusters and present some simple tests to show that 
these selection criteria are robust. 



3.1 Red Sequences 

With the unambiguous stars removed from our catalogs, 
we plot the colour-magnitude diagrams (CMD) of the re- 
maining objects and fit their red sequences. We present an 
overview of the red sequence fitting results here but for a 
more detailed description of the fitting procedure we refer 
the reader to Appendix [X] Briefly, we use the biweight esti- 
mator of Beers et al. (1990) to fit the colour- magnitude rela- 
tion (CMR), employing a statistical background subtraction 
method similar to that of Pimbblet et al. (2002). 

To check the quality of our photometry and red se- 
quence fitting procedure, we plot in Figure [3] the colour of 
the red sequence ((<?' — r') rs ) at a fixed M r i versus redshift. 
The ((?' — r') rs colour shown is calculated using a bootstrap 
method, finding the average (<?' — r') of the red sequence 
member galaxies lying within Am r / = ±0.5 mag on either 
side of the m r i value that separates dwarfs from giants (see 
Sj4] for detail of dwarf and giant selection). The strong de- 
pendence of (</ — r') rs colour on redshift seen in Figure [3] 
is a result of the 4000 angstrom break spectral feature, that 
is prominent in early type galaxies, shifting across the g' 
band with increasing redshift. The errors shown, which are 
typically smaller than the symbol size, account only for sta- 
tistical variance and do not include any systematics. We also 
show, as a dashed line, the expected (<?' — r') rs for our model 
red sequences (see appendix [XJ). The slope of the (</ — r') rs 
evolution is similar between the observations and the model, 
with a flattening at z ~ 0.4. There is however a significant 
systematic offset between the two from 0.05 < z < 0.3, 
as discussed in appendix [X] and is likely a consequence of 



the overly simplistic star formation history assumed in the 
model. 



3.2 Magnitude Limits 

We measure the ratio of luminous (giant) to faint (dwarf) 
red sequence galaxies (defined as \A(g' —r')\ < 0.2 from the 
best fit red sequence) using the absolute magnitude cuts de- 
fined by de Lucia et al. (2007) My = [-20.0, -18.2] which 
we transform to Megacam M r i using the single stellar pop- 
ulation red sequence model described in appendix [XJ Thus, 
we define giant galaxies as those with M r i < —20.32 and 
dwarf galaxies as those with —20.32 < M r i < —18.52. Note 
that these converted limits are identical to those used in 
Capozzi et al. (2010), who looked at the GDR using similar 
filters in the Sloan digital sky survey (SDSS). The excellent 
agreement with the Capozzi et al. (2010) limits at all red- 
shifts indicates that the difference introduced by comparing 
across photometric systems is negligible. 

At this stage, a choice in methodology must be made 
before these cuts can be applied to the data. One can either 
transform the data to absolute magnitudes or transform the 
absolute magnitude cuts to apparent magnitudes. We choose 
to do the latter because the former method requires a K- 
correction to all of the data, and since a if-correction based 
on only the g' and r' filters becomes degenerate beyond 
z ~ 0.6 this would allow additional high-redshift background 
objects to enter the dwarf and giant samples. Therefore our 
dwarf/giant sample selection is performed in observed space 
by transforming the above-mentioned absolute magnitude 
cuts following the relation: 

m r , = M r ,+51og(Di)+25-(K + E) (1) 

where Di represents the luminosity distance of the cluster 
in Mpc and the K-axiA E terms represent the if-correction 
and the evolution correction respectively. These corrections 
are described in detail in EI3.3I 

3.3 K + E Correction 

We assume that the luminosity evolution correction (E) 
takes the form of 2.5 x log(l + z), which is appropriate for 
passively evolving elliptical galaxies on the red sequence. 
The 7f -correction is calculated from the Chariot & Bruzual 
(CB07) synthetic spectra, assuming a solar metallicity Zq = 
0.02 with a 3% contribution by stellar mass from ultra-metal 
poor stars (see appendix |A"|) . The age of the spectral tem- 
plate is assigned by matching its redshifted (</ — r') colour 
to the (g' — r') colour of a galaxy on the best fitting red 
sequence at a given m r i . The appropriate K- correct ion is 
then calculated and this procedure is repeated for a sample 
of m r i values along the red sequence. Inverting equation [T] 
we generate the function that describes the if-correction as 
it depends on absolute magnitude. This is then interpolated 
for the dwarf and giant absolute magnitude limits described 
above to obtain the JsT-corrections that are specific to the 
colour-magnitude coordinates of said limits along the best 
fit red sequence. 

To verify that our Jf-correction is working properly 
we use a bootstrap method to compute the colour of the 
red sequence member galaxies in an m r i bin that extends 
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Figure 5. The giant-to-dwarf colour difference (GDCD) as it 
depends on redshift. Large symbols represent the error- weighted 
mean GDCD in each redshift bin. As a differential measurement, 
GDCD is unaffected by photometric zeropoint errors, making it 
a good benchmark for models of the red sequence. 

±0.5 around the cut in m r i that divides the dwarfs from 
the giants. Figure [4] shows the redshift dependence of this 
colour after if-correction. Once the Tf-correction has been 
accounted for the (g' — r') colour around the dwarf/giant di- 
vision remains relatively flat. There is one cluster (RXJ1524) 
that is a significant outlier, at z ~ 0.5 with a 7f -corrected 
(9' ~ r ')ra ~ 0.4. Inspection of the colour-magnitude dia- 
gram for this cluster reveals a second red sequence located 
at slightly bluer colour than that of RXJ1524, which is con- 
sistent with being the source of the observed colour bias. 
This contamination likely also biases the GDR measurement 
for this cluster but it does not affect our conclusions regard- 
ing the GDR of the ensemble. The constant colour seen in 
Figure [4] is expected for massive early type galaxies in this 
redshift range and lends confidence to the if-correction pro- 
cedure. The direct impact of the .ff -correction on GDR is 
investigated in t]4.3l 

3.4 Giant-to-Dwarf Colour Difference 

As a further check, we compute the (g' — r') colour of both 
dwarf and giant galaxies and examine the evolution of the 
giant-to-dwarf colour difference (GDCD) as defined by: 

GDCD( g i_ r i) = <( S ' - r') g ) - {(«,' - r') d ) (2) 

where ((</ — r') g ) and ((</ — r')d) are the mean colours of 
giant and dwarf galaxies respectively as calculated for indi- 
vidual clusters. Because the GDCD is a relative measure- 
ment, it is not affected by systematics that induce a shift in 
the photometric zero point over the entire image. Thus, it 
may provide a useful test of the galaxy properties in cluster 
formation simulations. Figure [5] shows the GDCD plotted 
against redshift. This quantity is effectively a measure of 



the red sequence slope as sampled by the individual mem- 
ber galaxies. From the model red sequences we expect the 
slope to increase with redshift out to z ~ 0.4 and then turn 
over as the (g' — r') filter combination becomes less sensi- 
tive to the 4000 angstrom break. Such behavior is reflected 
in the GDCD, though with a significant scatter. This high- 
lights the concern that the level of contamination by the blue 
cloud is a function of redshift and could falsely contribute to 
the evolution signal in the GDR. We attempt to correct for 
this effect in but we note here that if it were artificially 
creating a redshift trend in the GDR, then we would expect 
to see similar GDR values at both z = 0.1 and z = 0.5. This 
is not the case. 



4 GIANT-TO-DWARF RATIO 

In this section we discuss issues that deal directly with the 
GDR measurements. First, in section i]4.1l we describe the 
statistical background subtraction method and corrections 
for biasing effects introduced by uncertainty in (g' — r') 
colour. This is followed in t)4.2l by the results of GDR mea- 
surements for individual clusters. As mentioned in N3.2I we 
use magnitude limits defined by de Lucia et al. (2007) and 
convert to m r i(z). Our criterion for colour selection of dwarf 
and giant galaxies on the red sequence is based on the offset 
in (g' — r') from the best fit red sequence (\A(g' — r')\ < 0.2). 
We further require as a spatial criterion that the projected 
distance from the BCG is R < 750 kpc. Unless stated other- 
wise, the union of these three criteria define our dwarf and 
giant selection and are referred to as the default selection 
method. Within the literature however, there is no clear con- 
sensus on a set of criteria that defines the dwarf and giant 
populations, owing in part to varying data quality, instru- 
mentation and redshift range of observation. To facilitate a 
more direct literature comparison, in EI4.3I we examine the 
sensitivity of the GDR to variations in the selection crite- 
ria described above. Lastly, in H4.4I we compare our GDR 
values to previous results. 



4.1 Background Subtraction and Colour-Error 
Correction 

We apply the selection limits discussed above to our CMDs 
to find the number of dwarfs and giants in each cluster. The 
resulting numbers, however, are contaminated by interlop- 
ers. To correct for this contamination we apply a statistical 
subtraction method that uses the catalogs generated from 
the CFHTLS Deep field images (see £|2. 3|) . For each cluster, 
we take the m r i and (g' — r') limits for dwarfs and giants and 
apply these same limits to the CMD of the merged Deep field 
catalogs. We then correct the cluster counts in the following 
way: 

n c = n D — n b — (3) 
Ah 

where n c is the background-corrected number of cluster 
galaxies, n is the observed number of galaxies, rib is the 
number of galaxies in the background catalog and A and A b 
are the area of the target and background samples respec- 
tively. 
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Figure 6. The effect of background subtraction and colour-error 
correction expressed as a fraction of the total corrected counts for 
dwarfs (blue triangles) and giants (red squares) plotted against 
cluster redshift. The correction fraction becomes slightly nega- 
tive for the high redshift dwarfs when the colour-error correction 
begins to dominate over the background subtraction. 



An additional correction is made to account for the 
bias introduced by the uncertainty in the observed (<?' — r') 
colour. For instance, an intrinsically narrow red sequence 
will be broadened colourwise, which can result in galaxies 
being scattered outside of the colour selection limits, espe- 
cially at the highest redshifts. It is also possible for cluster 
members in the blue cloud to be observationally scattered 
into the red sequence, increasing the resulting number of red 
sequence galaxies and acting in opposition to the previously 
described effect. Both are statistically significant effects and 
failure to account for this observational bias can mimic an 
evolution in the GDR. 

To correct for the colour-error bias, we estimate the 
number of galaxies observed to be red sequence members as 
a fraction of the number that are intrinsically on the red se- 
quence (ie. the selection fraction). We simulate this process 
by taking the distribution of m r i and (g'—r 1 ) and perturbing 
each value by a normally distributed random number that is 
scaled to the la measurement uncertainty on each of these 
quantities. We compute the number of dwarfs and giants 
both before and after this perturbation and repeat the pro- 
cedure for 1000 realizations. We assume that the selection 
fraction (/ s ) can be approximated by f s ~ n rs /n' rs where 
n ra and n rs are the number of red sequence galaxies before 
and after the perturbation respectively. This approximation 
is valid as long as the difference (n rs — n' rs ) is small com- 
pared to n r3 . For the giants and dwarfs below z < 0.3 the 
effect is of order 5%, beyond z > 0.3 it grows to a maximum 
of 20%. The appropriate selection fraction is calculated for 
both the dwarf and giant populations and is compensated 
for by applying f s to n c from equation [3] to compute the 
total number of galaxies n t = n c / f s respectively. 



Figure 7. The Giant-to-Dwarf ratio (GDR) plotted against red- 
shift. Individual clusters are marked by the red open symbols. 
The blue line show the best fit linear redshift dependence. 

To assess the importance of background subtraction and 
the colour-error correction we show in Figure [6] the fraction 
of galaxies (/ c ) that must be subtracted to correct for these 
effects as a function of redshift. We define this parameter as 
follows: 

fc = (n - n t )/n t (4) 

The blue triangles and red squares in Figure [5] indicate 
dwarfs and giants respectively. The correction fractions are 
typically larger for the dwarfs than they are for the giants. 
The giants correction is typically constrained between 0-20% 
while the dwarfs exhibit a larger range in correction frac- 
tion. Note that f c can become negative at high z for some of 
the dwarf populations, which is expected due to increasingly 
fewer background galaxies in the selection box and increas- 
ingly larger colour-error; ultimately allowing the colour-error 
correction to dominate over the background subtraction. We 
show only the clusters with background and area-corrected 
number dwarfs plus giants larger than 50 (rid + n g > 50). 
Four clusters that are poorer than this limit can have f c val- 
ues in excess of 50% and we have excluded them from the 
rest of the analysis that follows. The excluded clusters are 
marked with an asterisk in Table [T] 

4.2 GDR results 

Figure [7] shows the GDR as a function of redshift for all 
of the clusters in our sample (filled red symbols). The data 
indicate that the GDR increases with redshift. The uncer- 
tainties include contributions from the shot noise in both 
the cluster field and background populations and a term 
that accounts for the error introduced by the presence of 
large-scale structure in the background. For each cluster we 
estimate the large-scale structure term by randomly placing 
100 apertures onto the CFHTLS Deep fields with radii of 
750 kpc and then calculating the 68% inclusion limits for 
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the number of background dwarfs and giants assuming a 
Poisson distribution. The solid blue line in Figure [7] shows 
the results of a best fit linear relationship between GDR 
and z parameterized by GDR = az + /3. We use the MP- 
FITEXY routine of Williams et al. (2010) to derive best fit 
values of a and /?, which uses the MPFIT package by Mark- 
wardt (2009). We also allow for some intrinsic scatter in the 
GDR {(Ti n t), which is allowed to vary in the fit. The data 
are bootstrap sampled for 1000 realizations to obtain an em- 
pirical estimate of the uncertainties on the parameters. We 
find a = 0.88 ± 0.15, /3 = 0.44 ± 0.03. The reduced x 2 of 
the fit is approximately equal to unity for all realizations, 
indicating that otnt is consistent with zero. If the large-scale 
structure error is neglected, however, the best fit parameters 
are similar a = 0.91 ±0.14, f3 = 0.47 ±0.03 but the required 
(Tint = 0.088±0.017. This suggests that large-scale structure 
is the dominant source of this scatter. A Spearman's rank 
test (non-parametric) of the GDR vs z correlation returns 
p — 0.54 with a significance of greater than 5a over the null 
hypothesis. The data and their best fit parameter values 
clearly favour an evolving GDR, but we reserve further dis- 
cussion of this and a detailed comparison with results from 
the literature for i]4.4l A list of the GDR measurements for 
each cluster, along with the number of dwarfs and giants 
both before and after the background±colour-error correc- 
tion (see M4 . 1 p is given in Tabled 

4.3 Sensitivity Analysis 

One of the issues that complicates the comparison of GDR 
(and dwarf-to-giant ratio) between studies is the varying 
selection criteria used to define the dwarf and/or giant pop- 
ulations. The _ff-correction method, the projected cluster- 
centric distance within which to include galaxies in the 
dwarf and giant populations as well as the maximum colour- 
residual with respect to the CMR that defines a red sequence 
member are examples of parameters that often vary between 
different studies. In this section, we compute the GDR using 
a range of selection criteria to see how much of the apparent 
discrepancy can be explained through these choices. 

To examine how our results depend on the adopted K- 
correction described above, we compare the GDRs obtained 
using two different i\-correction methods. The results are 
plotted in Figure [HJi. The first Tf-correction method (K0) 
is the one outlined above using the CB07 models and as- 
sumes solar metallicity and a template age that is deter- 
mined by matching the observed (</' — r') colour. The next 
method (K2) uses the Tf-correction code of Chilingarian, 
Melchior and Zolotukhin (2010) which has been shown to 
agree with the empirically derived corrections of Blanton 
& Roweis (2007) up to z=0.5. We also show the results of 
neglecting the ^-correction completely (no K). 

The differences in GDR between one _R"-correction 
method and another are typically smaller than the GDR 
error in an individual redshift bin. Moreover, the selection 
of either K-correction method does not affect the conclusions 
regarding the presence or absence of a trend with redshift, 
only how strong this evolution might be. However, neglect- 
ing the if-correction would lead to a contradictory view of 
the redshift evolution in GDR. 

Although it is not shown in Figure [HJi, we also test the 
effect of neglecting the evolution correction (ie. E = 0, see 



i]3.3|l . We find that contrary to neglecting the A'-correction, 
this leads to a significant steepening of the GDR vs z rela- 
tion (a = 1.61 ± 0.15). Therefore, to obtain a flat GDR vs z 
by varying the _E-correction would require E to be param- 
eterized with a significantly stronger redshift dependence 
than what we have assumed here. Such corrections are un- 
reasonable for the early-type galaxies that make up the red 
sequence. 

Figure [Hp shows the effect on GDR of varying the aper- 
ture size (rum), centred on the BCG, within which to count 
red sequence galaxies. The data shown in the figure are for 
galaxies within 500 kpc, 750 kpc, 1000 kpc, 1250 kpc and 
1500 kpc. There is some variation in the results, most no- 
tably in the lowest redshift bin, but it is all within the error 
in an individual redshift bin and thus the measurements 
based on the various apertures are in agreement. We note 
that the largest apertures show the smallest GDR errors. 
This is expected since they are based on Poisson statistics 
and the GDR is a ratio. As the number of galaxies in the 
population increases, the error as a percentage on that num- 
ber decreases. We argue, however, that it is better to use a 
smaller aperture such as ru m = 750kpc because the po- 
tential for systematic error introduced by the background 
subtraction procedure, which is not included in the error 
estimates shown in Figure [Hp, grows as the aperture size in- 
creases. It is also interesting to examine the differences in 
the GDR between the cluster centre and cluster outskirts. 
Such a discussion does not fit well in this section however, 
which is intended to assess the possible variation in GDR 
introduced by commonly used selection criteria. For a dis- 
cussion of the GDR in the inner vs outer parts of clusters 
we refer the reader to ij6] 

We consider the possibility that BCG position may be a 
poor indicator of the cluster centre as a function of redshift. 
Depending on the radial profile of the GDR, this effect could 
potentially mimic GDR evolution. To test this, we measure 
the GDR in 0.05 < z < 0.15 clusters using a centre position 
that is offset between to 6.5 arcminutes from the BCG. 
This shift corresponds roughly to 750 kpc at z = 0.1. The 
results are consistent with those obtained using the BCG 
as the cluster centre, indicating that the GDR evolution is 
unaffected by centreing errors on these angular scales. 

We show in Figure [HJ: how varying the colour selection 
criteria affects the measurement of the GDR. Here we com- 
pare several cases of two distinct selection criteria. The first 
requires galaxies to be within \A(g' — r')\ < 0.2, which is 
the assumed criteria throughout the rest of this work, but 
we also show here the results for \A(g' — r')\ < 0.1 and 
\A(g' — r')\ < 0.3. The second selection method that we in- 
vestigate is based on the observed scatter in (g' —r') of the 
red sequence (<r), assuming a gaussian distribution about 
the best fit colour-magnitude relation. Figure [HJ: shows the 
GDR of galaxies that lie within \A(g' — r')\/a < 3 and 
\A(g' — r')|/<7 < 2. The effect of colour selection criteria 
is larger than that introduced by the ru m selection crite- 
ria and comparable to the effect of varying the Tf-correction 
method. However, the change in GDR is always smaller than 
the error in an individual redshift bin, which demonstrates 
that the choice of a particular colour criterion does not affect 
conclusions about the presence of evolution. 
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Table 2. The GDR measurements. The columns give the name, the redshift (z), the fully-corrected GDR (GDR), the statistical 
uncertainty on the GDR (ctgdr), the number of giants and dwarfs after data corrections (n g & respectively), the number of giants 
and dwarfs prior to data corrections (n g aw & n ^ aw respectively), the maximum m r i limits used for defining the giant and dwarf samples 
(rn 9 r , &i mr, respectively), the colour of the red sequence ((g' — r') rs ± 0.01) - see i|3.1H and the X-ray temperature (T x ) as measured 
with ASCA (Horner 2001). Additional temperature measurements are added from BAX and marked with a * in the last column. The 
numbers of galaxies have been adjusted to account for the image area that was masked in the detection procedure. 
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4.4 Comparison with Literature 

The disagreements in the literature regarding the evolution 
of the GDR/DGR are primarily driven by the interpreta- 
tion of various results and not necessarily by disagreements 
among the data themselves. In a study based on DGR data 
culled from the available literature, Gilbank and Balogh 
(2008) show that if one neglects the z > 1 clusters of An- 
dreon (2008), where the filters do not bracket the 4000 



angstrom break, DGR vs z data is fit reasonably well by 
a single power law in (1 + z). Andreon (2008) puts forth an 
important caution that one must be careful to account for 
the presence of an intrinsic scatter in DGR from one clus- 
ter to another so as not to be biased by outliers. Andreon 
(2008) reports an intrinsic scatter <Ji nt = 0.13 ± 0.06, ar- 
guing that this is enough to reject trends as steep a ~ 1.3 
found in de Lucia et al. (2007), where ai nt was not taken 
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into account. The author also provides a caveat that the 
observed spiral morphologies of some nearby red sequence 
galaxies (eg. Butcher & Oemler 1984) indicate recent star 
formation, colour transformation and red sequence buildup. 
However, the associated GDR evolution, if any, is not steep 
enough to be detected in small samples of ~ 25 clusters. To 
contribute to this discussion we compare these result with 
our own for ~ 100 clusters. 

As discussed in g^l we find a = 0.88 ± 0.15, p = 
0.44 ± 0.03 and Oi n t consistent with zero. If the large scale 



structure term is neglected in the GDR uncertainty we ob- 
tain a ai n t = 0.088 ± 0.017 which is consistent with the 
value inferred by Andreon (2008). Our best fit slope is sig- 
nificant, though not as steep as that found in de Lucia et 
al. (2007) . Figure [5] shows how our results compare with 
those from the literature. The large square symbols indicate 
the error- weighted mean GDR for our sample in several red- 
shift bins. The solid blue line shows our best fit linear re- 
lation described above. Data taken from the literature are 
overplotted. The Gilbank et al. (2008) results are based on 
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Figure 8. (a) The giant-to-dwarf ratio using various methods of X-correction. The different symbols correspond to 2 different K- 
correction methods (see text for details) as well the effect of neglecting the ^-correction, (b) The giant-to-dwarf ratio using various 
apertures sizes for red sequence membership. The different symbols refer to galaxies within r ;im in projected distance from the BCG 
with ri im =500 kpc, 750 kpc, fOOO kpc, 1 250 kpc and 1500 kpc. (c) The giant-to-dwarf ratio as computed using various colour selection 
criteria. The different symbols show \A(g' - r')\ < 0.1, \ A(g' - r')| < 0.2, A(g' - r')| < 0.3, \A(g' - r')\ja < 2 and \ A(g' - r')\/a < 3. 
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A de Lucia et al. 2007 
n Andreon 2008 
& Gilbank et al. 2008 (corrected) 
o Capozzi etal. 2010 




Figure 9. The ratio of the number density of giant galaxies to 
that of dwarf galaxies as a function of rcdshift. The data has 
been binned in redshift bins of Az = 0.05 (red squares). We also 
show various data from the literature. Note that the Gilbank et 
al. (2008) values have been adjusted to be comparable directly 
with ours (see text for description). 



slightly different magnitude limits for defining dwarfs and 
giants than those used here. To make a fairer comparison 
we compute the GDR by applying their magnitude cuts to 
our data, then calculate a correction based on the difference 
between the best fit a and /3 from those obtained using our 
default cuts. We apply this correction to the Gilbank et al. 
(2008) values and plot the 'corrected' values in Figure [9] 
There is very good agreement between our results and the 



'corrected' values of Gilbank et al. (2008) despite the differ- 
ence in redshift range and cluster masses probed. Our GDR 
values are also in agreement with the low redshift GDR re- 
sults of Capozzi et al. (2010), which are based on SDSS clus- 
ter catalogs and the de Lucia et al. (2007) magnitude limits. 
When compared directly to the de Lucia et al. (2007) results 
our data do not seem to agree. However, there is an impor- 
tant difference between the two analyses. The magnitudes 
used in this paper are taken from SExtractor's MAG_AUTO 
parameter, which attempts to follow the shape of the galaxy 
surface brightness profiles. The de Lucia et al. (2007) mag- 
nitudes are extracted using apertures of 1.0", significantly 
smaller than the sizes within which our magnitudes are cal- 
culated. Using a 1.0" fixed aperture truncates the galaxies 
at smaller radii, effectively reducing their fluxes. This im- 
parts a stronger bias against the giants than dwarfs, since 
the former are intrinsically more extended on the sky lead- 
ing to a reduction in the value of the GDR and a steepening 
of the GDR vs z relation. We are not able to fully reproduce 
an analogous GDR measurement to that of de Lucia et al. 
(2007) because of the different angular resolution scales of 
CFHT and HST. As an instructional exercise we remeasure 
our GDR values using a 3.0" aperture size for galaxy magni- 
tudes, which is more suitable given the spatial resolution of 
CFHT. We verify that this leads to a significant reduction in 
the GDR values compared to our default extraction method. 
The resulting best fit to the fixed-aperture GDR vs z rela- 
tion has a slope a — 1.90, which represents a steepening of 
Aa ~ 1. We also note that because this effect operates more 
noticeably at lower redshifts as the angular sizes of galaxies 
grow and the truncation becomes more extreme, it may pro- 
vide an explanation for the very steep redshift trend seen by 
de Lucia et al. (2007). 
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Figure 10. GDR versus 2 with clusters coded by X-ray temper- 
ature information. The solid black line shows the best fit linear 
relation found in i|4.2l The GDR evolution can be represented by 
a single function of redshift regardless of T x for the full range of 
temperature and redshift probed by our sample. 
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Figure 11. The number of giants (red triangles) and dwarfs (blue 
squares) as a function of cluster T x for clusters in the range 0.1 < 
Z < 0.3. Solid and dashed lines correspond to the best fit relations 
for dwarfs and giants respectively. 



5 DEPENDENCE ON CLUSTER MASS 

As we discuss in the previous section, the GDR data can 
be modelled by a linear trend with z. After accounting for 
large-scale structure in the background, we find no intrinsic 
scatter. This is qualitatively consistent with the results of 
Lu et al. (2009) who find that the large-scale structure un- 
certainty comprises a significant portion of the uncertainty 
in DGR for individual clusters. This result is consistent with 
a GDR that is a function of redshift only and is indepen- 
dent of any other global cluster properties, at least over the 
range of redshift and cluster mass investigated in this paper. 
Several authors in the literature, however, have presented 
evidence for a dependance of GDR on global properties re- 
lated to cluster mass. For instance, de Lucia et al. (2007) 
find a higher GDR in their a > 600 km/s sample than that 
in their a < 600 km/s sample, with a difference in GDR of 
about 0.1 to 0.2. They also analyze data from the SDSS for 
clusters at z ~ and find that this a dependence is reversed. 
Gilbank et al. (2008) find an elevated and flatter GDR for 
poorer clusters. Capozzi et al. (2010) argue that there may 
be a weak correlation between the GDR and cluster X-ray 
luminosity. Furthermore, as the GDR is different in the field 
(Gilbank & Balogh 2008), some mass dependence might be 
expected. Because our cluster sample is not a representative 
sample, it may be that some portion of the observed GDR 
evolution arises as a consequence of potential covariance be- 
tween redshift and cluster mass (more high mass systems 
at high z). However, given the good agreement between our 
GDR results and those based on cluster samples with lower 
mean mass (eg. Gilbank et al. 2008 and Capozzi et al. 2010) 
we do not expect this to be the case. Nonetheless, in this 
section we investigate this possibility in more detail. 

Many of the previous results that indicate a dependence 



on cluster mass are based on using cluster richness as a mass 
proxy. One problem with using richness is that there is a 
large scatter in the mass-richness relation of up to ±150% 
in mass at a given richness (Rozo et al. 2009). Moreover, the 
richness, which is approximately the number of dwarfs plus 
giants (rid + n g ) and the GDR (n g /nd) are by construction 
naturally covariant. A better proxy for cluster mass is the X- 
ray temperature of the hot intra-cluster medium (T x ), which 
predicts the mass with significantly lower scatter (36% in 
mass at a fixed T x inside rsoo) and is independent of galaxy 
counting. 

A subset of the clusters in our sample have been ob- 
served with the Advanced Satellite for Cosmology and As- 
trophysics (ASCA) and analyzed by Horner (2001). For 
maximum consistency we take T x values from the Horner 
(2001) catalog for as many of our clusters as possible (59 
systems). We supplement this with an additional 15 T x mea- 
surements from the BAX X-ray galaxy cluster databas^fl 
giving us a total of 72. A complete listing of the T x data is 
given in Tabled 

To investigate the effect of T x on GDR evolution we 
show in Figure [10] the plot of GDR vs z after separating 
the data into 3 different T x categories. The data points in 
Figure [10] coded by T x > 6keV (red circles), T x < 6keV 
(blue snowflakes) and no T x information (black crosses). Er- 
ror bars have been omitted for clarity but are listed in Table 
[2] Figure [10] shows that indeed there are more high T x sys- 
tems at the high z end of our sample and conversely more 
low T x systems at the low z end. Despite this selection ef- 
fect however, the data in each T x category are consistent 
with the same best fit GDR vs z relation found in section 
14.21 Specifically, we find that the best fit parameters are 

3 http://bax.ast.obs-mip.fr/ 
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(a = 0.90 ± 0.31, P = 0.49 ± 0.07) for the high T x clusters, 
(a = 0.92±0.27, (3 = 0.51 ±0.05) for the low T x clusters and 
(a = 0.98 ± 0.22, /3 = 0.40 ± 0.04) for the clusters with no 
T x information. This consistency is compatible with a GDR, 
and GDR evolution that do not depend on cluster mass over 
the range of mass and redshift probed by our sample. 

A subset of 27 clusters in our sample are part of a more 
detailed campaign to measure the masses using a combined 
X-ray + weak lensing method. We compare the GDR di- 
rectly to the weak-lensing masses (taken from Hoekstra et 
al. in prep.) but find no evidence for a correlation. This is 
consistent with the result obtained using T x as a mass proxy. 
From the X-ray portion of the analysis we obtain central en- 
tropy values (So) for the subset (taken from Mahdavi et al. in 
prep.). The value of So is often used to discriminate between 
cool-core and non-cool-core clusters, so it is interesting to 
examine if/how the GDR varies with So- We find no cor- 
relation between these two parameters, however, indicating 
that the GDR is similar between cool-core and non-cool-core 
clusters. Lastly, we consider the GDR as a function of the 
projected distance between the BCG and the peak of the 
cluster X-ray emission, an indicator of the dynamical state 
of the cluster. Using values from Bildfell et al. (2008) we 
find no evidence for a correlation between GDR and BCG 
to X-ray peak distance. These results are all consistent with, 
and expected for a GDR vs z relation that shows no intrinsic 
scatter. 



6 DRIVER OF EVOLUTION 

The physical mechanism(s) responsible for the evolution in 
the GDR for 0.05 < z < 0.55 is not evident. The two main 
processes affecting red sequence galaxies often discussed in 
the literature are dry mergers, which reduces the number of 
galaxies on the red sequence and late-time quenching of star 
formation, which increases the number of galaxies on the 
red sequence. Both of these mechanisms are taking place in 
and around clusters but the extent to which they operate 
and how their behavior may change with redshift remains 
elusive. Bell et al. (2004) and Faber et al. (2007) argue that 
both quenching and dry merging are required to explain 
observations of the red sequence buildup since z ~ 1. Ob- 
servational evidence for dry mergers of early type galaxies 
at < z < 0.8 is presented in van Dokkum (2005) and Tran 
et al. (2005). On the other hand, work by Cimatti et al. 
(2006) and Scarlata et al. (2007) argue that the buildup of 
the red sequence since z ~ 0.7 and z ~ 1, respectively, can 
be explained by the late-time quenching of star formation 
alone. With the exception of Tran et al. (2005) however, 
these studies are based on samples containing galaxies in 
the field and not exclusively on cluster galaxies where the 
high-density local environment may be important. To inves- 
tigate this further, we examine the evolution of the dwarf 
and giant populations separately. 

We showed in [J5]that the GDR is not strongly depen- 
dent on cluster mass (or T x ). When considered indepen- 
dently however, n g and rid are both expected to scale with 
cluster mass; similar to the known scaling between richness 
and T x (Yee & Ellingson 2003). Such a scaling is implied by 
the existence of a mass-richness relation. This effect must be 
accounted for before we examine the evolution of the dwarf 



and giant populations separately. We show in Figure [TT] the 
dependence of galaxy number n g and rid on T x . The data are 
fit using a linear relation between number density and T x . 
To minimize the effect of the correlation between T x and z 
imposed by our sample selection, we only fit clusters in the 
range 0.1 < z < 0.3. We find a best fit slope a g = 6.9 keV -1 
for the giants and ad = 9.8 keV -1 for the dwarfs. Based on 
this temperature dependence we correct all n g and rid values 
to a common T x = 6 keV. 

We plot the corrected values as a function of redshift in 
Figure 112b ,. The data are combined using a constant num- 
ber of clusters in each bin (horizontal error bars indicate bin 
limits). Both the number of dwarfs and the number of giants 
increase from z ~ 0.55 up to z ~ 0.3, with the dwarfs show- 
ing the more dramatic evolution. From z ~ 0.3 to z ~ 0.2 the 
numbers of dwarfs and giants remain relatively constant. Be- 
low z ~ 0.2 the number of dwarfs remains constant but the 
number of giants begins to decrease significantly. The results 
suggest that the GDR evolution at z < 0.2 is dominated by 
mergers between giant galaxies. When considered with GDR 
results at higher redshift, the data support a transition in 
the dominant mechanism governing GDR evolution, which 
is often attributed to quenching in dwarfs. This interpreta- 
tion could also explain the findings of Lu et al. (2009) , who 
measure a rapid DGR evolution at < z < 0.2 and then 
a flattening of the relation out to z ~ 0.4. Though not dis- 
cussed in their work, Figure 19 of Lu et al. (2009) clearly 
shows a reduction in the number of giants from z ~ 0.2 to 
z ~ 0.1. 

Given the evidence supporting merger-driven evolution 
at low redshift, it is natural to ask whether this result is com- 
patible with theoretical predictions. The orbit of a galaxy 
about the cluster centre will decay with time due to dynam- 
ical friction. The relevant timescale for the orbital decay of 
an individual galaxy can be expressed as: 



tdyn 



6 x 10 yr 



lOOOkm/s I V 0.25Mpc 



(M/L V )gal 



WMq/Lq 



(5) 



with line-of-sight velocity dispersion ay, cluster core radius 
r c , galaxy mass m and the galaxy's V-band mass-to-light 
ratio (M/Ly). For a giant galaxy with m = m* and a 
M/Ly — WMq/Lq orbiting near the centre of a cluster 
where oy ~ 600 km/s and r c = 0.25 Mpc the expected 
timescale is td yn ~ 3.6 Gyr. Over the redshift range of our 
data (0.05 < z < 0.55) the total change in cosmic time 
is At age = 4.7 Gyr, which we use as an estimate of the 
elapsed interaction time. This is considered to be a conser- 
vative estimate because these clusters are likely assembled at 
much earlier epochs. Since td yn < At age we should indeed 
expect to see the signature of the most massive galaxies 
merging together on these timescales. This prediction is in 
good qualitative agreement with the decrease in the number 
of giants galaxies at low redshift in Figure 112b . Further- 
more, quenching is expected to become less important at 
late times due to the declining specific star formation rate 
(at fixed galaxy mass) in newly accreted field galaxies (eg. 
Juneau et al. 2005). 

If the GDR evolution seen at low redshift is driven by 
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Figure 12. The evolution of the number of dwarfs (blue squares) and giants (red triangles) after correcting to a common T x = 6 keV. 
Vertical error bars show the uncertainty on the error-weighted bin centroid and horizontal error bars indicate the width of the bins. The 
data in the left panel (a) shows the number of galaxies in a projected radius r < 750 kpc (default region) while right panel (b) shows 
the number of galaxies in a region 750 kpc < r < 1500 kpc. 



orbital decay and mergers in the giant population then we 
might expect to see less of this effect taking place in the 
cluster outskirts, where the density of the ambient medium 
is lower. To test this we plot in Figure 112b the number of 
giants and dwarfs (adjusted to T x = 6 keV) in an annulus 
with the range of projected radii 750 kpc < r < 1500 kpc. 
The number of giants rises from z ~ 0.5 and peaks near 
z ~ 0.25 declining and flattening below z < 0.15. The num- 
ber of dwarfs generally increases from z ~ 0.5 to z ~ 0.05 
but the distribution is somewhat noisier than similar data at 
smaller projected cluster-centric distance (Figure Il2k). The 
triggering of star formation as galaxies fall into the cluster, 
pre-processing of galaxies in groups, and the increased con- 
tamination by interloper galaxies may all affect the number 
counts in the cluster outskirts and complicate their inter- 
pretation considerably. 

It is difficult to explain all of the features of Figure [12] 
in the context of mergers and quenching alone. The rate 
of newly infalling cluster galaxies, as it depends on redshift 
and galaxy mass, is also likely to be a key ingredient. The 
merger tree models of de Lucia et al. (2011) indicate that 
up to ~ 50% of the massive galaxies (log M*/M® > 10) 
that fall into clusters do so after z = 1. Taking this into 
consideration, it may be more appropriate to describe the 
GDR evolution at low redshift as competition between the 
galaxy infaU rate and the effects of dynamical friction, with 
quenching as a secondary effect. 

For a final check we compare in Figure [T3] the GDR 
evolution in the cluster outskirts (750 kpc < r < 1500 kpc) 
with that in the cluster centre (r < 750 kpc). The data 
indicate that the GDR evolution, as measured in the clus- 
ter interior and outskirts, are in good agreement with each 
other. This suggests that the physical mechanism(s) respon- 
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Figure 13. The evolution of the GDR in two distinct regions 
of projected-radius. Filled symbols represent the GDR inside a 
projected radius limit r < 750 kpc (default region), while open 
symbols represent the GDR in the range 750 kpc < r < 1500 kpc. 



sible for driving the GDR evolution operate both at large 
and small cluster-centric distance. This may simply reflect 
the fact that clusters are assembled from groups and that 
infalling galaxies are likely to have already undergone some 
processing in high-density environments. 
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6.1 Integrated Stellar Masses 

To test the interpretation of a merger-driven GDR evolution 
at z < 0.2, we measure the integrated stellar mass of giants 
(Mi) and dwarfs (M«), separately, as a function of z. We 
rearrange equation[T]to obtain (_K" + 15)-corrected M r i values 
for individual galaxies which are then converted to r'-band 
luminosities L r assuming an r'-band absolute magnitude for 
the sun M ri @ = 4.650- The luminosities of interloper galax- 
ies are calculated in a similar way using the CFHTLS Deep 
background populations described in £14- 1 1 The integrated 
dwarf or giant luminosity is obtained for each cluster by 
summing L r values for a given population and then subtract- 
ing the appropriate background. We also apply the colour 
error correction described in H4.ll To convert the luminos- 
ity to a stellar mass we assume a universal stellar mass-to- 
light ratio in the r'-band of M/L r = 1.48, which is obtained 
from the scaling relation of Zibetti et al. (2009) for elliptical 
galaxies with (g — r) = 0.61. We note however, that the ex- 
act normalization is unimportant for our purpose, as we are 
primarily interested in the relative change in stellar mass. 

To make a fairer comparison of the masses, we fit the 
log M* versus T x relation (for 0.1 < z < 0.3) and correct to a 
common cluster temperature of T x = 6 keV. This procedure 
is analogous to that in ijH regarding the number densities. 
We show in Figure [14] the T^-corrected values of Ml and 
Mf as a function of redshift. The data reveal that the mean 
stellar mass in dwarf galaxies ((M,)) increases from z ~ 0.55 
to z ~ 0.25 and remains roughly constant for z < 0.25. In 
contrast, the mean stellar mass in giants ((A/*)) is roughly 
constant from z ~ 0.55 to z ~ 0.2 and then declines at 
z < 0.2. 

In [JSjwe propose that the decrease in the number of gi- 
ants at z < 0.2 can be explained by a significant number of 
giant-giant mergers (merger scenario). If this is the case then 
one might expect (M») to remain constant over this redshift 
range. However, there are several reasons why mergers may 
act to deplete the integrated stellar mass. Tidal features cre- 
ated in mergers are believed to contribute to the production 
of intra-cluster stars (ICS), which are not gravitationally 
bound to any particular galaxy and make up ~ 10% — 30% 
of the stellar mass in clusters (Gal- Yam et al. 2003; Gonza- 
lez et al. 2005; Zibetti et al. 2005; Krick & Berstein 2007, 
Sand et al. 2011). Furthermore, the existence of a tight rela- 
tionship between effective radius and mean effective surface 
brightness (Kormendy 1977) implies that as elliptical galax- 
ies grow, they become increasingly diffuse. Consequently, the 
fraction of a galaxy's stellar mass that is below our detection 
threshold may be a function of the number of mergers that 
galaxy has undergone. 

We see from Figure [T^] that there is a 29±13% reduction 
in the number density of giants from z ~ 0.2 to z ~ 0.05. 
The corresponding reduction in (Ml) over this period is 
38±14%. Assuming the merger scenario is correct, a sizeable 
fraction of the stellar mass involved in a giant-giant merg- 
ers, up to 100% in some cases, must be shifted to the ICS 
or drop below our detection threshold. As a caveat to these 
statistics, Gonzalez et al. (2007) finds that the ICS fraction 
scales inversely with the velocity dispersion of the cluster. 
If ICS production is much more efficient in the lowest mass 

4 http://mips.as.arizona.cdu/ cnaw/sun.html 
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Figure 14. Integrated masses of galaxies associated with the 
dwarf population (blue squares) and giant population (red trian- 
gles) as a function of redshift. The data have been scaled to a 
common cluster Xray temperature T x = 6 keV. Individual clus- 
ters are shown as smaller, lighter symbols, while binned data is 
shown as larger, darker symbols. Error bars on binned data indi- 
cate the bin width (horizontal) and the standard deviation within 
the bin (vertical). The mass in dwarf galaxies is roughly constant, 
while the mass in giants is decreasing below z < 0.2. 

clusters then this effect may be exacerbating the observed 
drop in (Mi) at z < 0.2, where our sample is dominated by 
clusters with low T x . To test this we repeat our measure- 
ment using only those clusters with T x > 6 keV. With this 
selection, the drop in (Ml) at z < 0.2 is less pronounced and 
(Ml) is consistent with a constant value over all z. These 
results are suggestive but their interpretation is limited by 
the large standard deviation within each bin. 

Because of the above mentioned considerations we are 
unable to confirm/deny the merger scenario with the inte- 
grated stellar masses alone. To make progress requires mea- 
surement of the ICS fraction for all of the clusters in our 
sample and better constraints on the redshift-dependance 
of the ICS fraction, both of which are outside the scope of 
this paper. 



7 CONCLUSIONS 

We analyzed a large sample of 97 galaxy clusters observed in 
g' and r' with CFHT-MegaCam. After carefully accounting 
for several observational biases and statistically removing 
interloper galaxies using g' and r' galaxy catalogs generated 
from the CFHTLS Deep fields as a control, we measured the 
GDR in galaxy clusters over the range 0.05 < z < 0.55. We 
find several interesting results, the most important of which 
we summarize below. 

We detect evolution in the GDR over the redshift range 
0.05 < z < 0.55, which can be parametrized by a simple lin- 
ear relation with z (GDR — az + j3). The best fit parameters 
are a = 0.88 ±0.15 and fj = 0.44±0.03 with an estimate for 
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the intrinsic scatter in the relation that is consistent with 
zero. Neglecting the uncertainty introduced by the presence 
of large-scale structure leads to a perceived intrinsic scatter 
of Ui n t = 0.088 ± 0.017, which is consistent with previous 
estimates (Andreon et al. 2008). 

Varying the method of selecting dwarf and giant galax- 
ies by altering the BCG-centric distance threshold, colour se- 
lection criteria or K- correct ion prescription does not change 
the GDR significantly. So long as one makes reasonable 
choices about the selection criteria and if-correction, then 
varying these parameters cannot hide the signature of an 
evolving GDR. 

The GDR values found here agree well with the those 
from the literature that are derived directly, solely from mea- 
surements of the number of luminous and faint galaxies. 
However, our GDR results are not in agreement with those 
that are inferred from parametric fits to the luminosity func- 
tion (eg. Andreon 2008, Crawford et al. 2009). 

The GDR does not correlate with cluster mass estimates 
via X-ray temperature or weak-lensing. Similarly, we find 
no evidence for a correlation between GDR and the central 
entropy of the ICM or the dynamical state of the cluster 
(BCG to X-ray peak offset). This is further validated by the 
good agreement between our GDR values and those found in 
studies based on lower mass clusters (eg. Gilbank et al. 2008, 
Capozzi et al. 2010). These result agree with our finding that 
the intrinsic scatter in GDR vs z is consistent with zero. 
When considered individually however, both the number of 
dwarfs and the number of giants are strongly correlated with 
T x , which is expected from the known correlation between 
richness and T x (Yee & Ellingson 2003). 

Separate inspection of the evolving number of dwarfs 
and giants, after correcting to a common T x , suggests a 
change in the primary physical mechanism responsible for 
the evolution in the GDR over the redshift range 0.05 < 
2 < 0.55. From z ~ 0.55 to z ~ 0.25 we observe an over- 
all increase in the number of both dwarfs and giants, with 
a more rapid increase among the dwarfs. Below z ~ 0.2 
however, the number of dwarfs is roughly constant, while 
the number of giants decreases significantly. To explain the 
transition at z ~ 0.2, we argue for a significant number of 
mergers within the giant population, caused by orbital decay 
due to dynamical friction. This is supported by calculations 
of the expected orbital decay time for giant galaxies, which 
yield timescales smaller than the span in look-back time of 
our sample. In this scenario, the GDR is governed by the 
competition between mergers, which become dominant at 
z < 0.2, versus the processes that act to increase the rel- 
ative number of dwarf galaxies on the red sequence, such 
as the late-time quenching of star formation and the red- 
shift dependence of the accreted-galaxy mass function. We 
consider the latter to be particularly important given the 
results of de Lucia et al. (2011) indicating that up to 50% 
of log M,/Mq > 10 galaxies fall into clusters at z < 1.0. A 
GDR evolution that is driven by a reduction in the num- 
ber of giants at z < 0.2 could also explain why Crawford 
et al. (2009) find no evolution in the faint-end slope of the 
luminosity function at low redshift. 

To test the merger-driven GDR scenario, we measure 
the integrated stellar mass in giants and dwarfs as a func- 
tion of redshift. We find that (M?) increases from z ~ 0.55 
to z ~ 0.25 and remains roughly constant for z < 0.25, while 



(Mi) is roughly constant from z ~ 0.55 to z ~ 0.2 and then 
declines at z < 0.2. Because of the large scatter present in 
the values of Ml at fixed z and the uncertainty regarding 
the production of ICS, it is not possible to confirm nor deny 
the merger scenario with this measurement alone. To make 
progress requires measurement of the ICS fractions of clus- 
ters in our sample and better constraints on the evolution 
of the ICS. 

The results presented in this paper clearly indicate an 
evolving GDR in our sample of 97 galaxy clusters at in- 
termediate redshift. We are cautious, however, not to over- 
interpret the data in light of the fact that the sample is 
incomplete. Correlation between the GDR and global clus- 
ter properties, along with a covariance between said prop- 
erties and z in our sample, could be misinterpreted as an 
enhancement of the GDR evolution. The data presented in 
ij5] however, along with the good agreement with the litera- 
ture values over a range of sample selection properties, and 
an intrinsic scatter that is consistent with zero, all indicate 
that this effect is negligible. A future study of the GDR using 
a representative sample of at least as many clusters, along 
with spectroscopic redshifts for identifying cluster members 
would be useful to confirm/deny these findings and allow 
for a clearer interpretation that is unaffected by potential 
selection biases. Numerical models of cluster formation that 
focus on the competition between mergers of existing clus- 
ter members versus the recent infall of new cluster galaxies, 
as it depends on galaxy mass and redshift, would be valu- 
able for assessing the proposed mechanisms for driving GDR 
evolution. 
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APPENDIX A: RED SEQUENCE FITTING 



As mentioned in i|3.1l we use the biweight estimator of Beers 
et al. (1990) to fit the colour-magnitude relation (CMR) for 
the galaxy distribution inside a projected radius of 500 kpc 
from the BCG. Lacking spectroscopic redshifts for the clus- 
ter galaxies, we begin the fitting procedure by estimating a 
region in the CMD that is likely to be minimally contami- 
nated by interloper galaxies. To the points lying inside this 
region we then apply a red sequence fitting procedure similar 
to that used by Pimbblet et al. (2002). We first determine 
the probability that a galaxy belongs to the field population 
using equation IA1I 



Pf(m r ',g' - r') 



T, f (m r ,,g' - r') 



Z f (m r ,,g' -r') + £ c (rrv,<?' -r') K ' 

where £/ is the surface density of field galaxies, E c is the 
surface density of cluster galaxies and both are functions of 
m r i and (<?' — r'). The surface densities are estimated by 
generating a 2-d histogram of galaxies in colour-magnitude 
space using a bin height and width of A(g' ~ r') = 0.078 and 
Am r i = 0.56. The denominator in equation lAll is calculated 
using galaxies within a projected radius of 500 kpc from the 
BCG where there is a mix of the cluster and field popula- 
tions, while the numerator is calculated using galaxies that 
are outside a projected radius of 1.5 Mpc where the sample 
should be more dominated by the field population. We chose 
the 1.5 Mpc limit in order to be as far away from the BCG 
as possible while maintaining a large enough field-sampling 
region at the low redshift end of the survey z ~ 0.05. At this 
redshift roughly 45% of the image is available for estimating 
the field population. With each galaxy assigned a probability 
of being a field member (Pfieid) and a probability of being 
a cluster member (1 — Pfieid) we then fit the red sequence 
to the points inside the selection box. We randomly reject 
points inside the initial selection box based on their prob- 
ability of being members of the field and then perform an 
unweighted, linear least-squares fit on the remaining data. 
This procedure is repeated 100 times. We take the median 
values of these 100 solutions as the best fitting slope and 
intercept parameters with their errors determined by the 
standard deviation. 

The results of the red sequence fitting for a subset of 20 
randomly selected clusters are plotted in Figure IAT1 In each 
panel we show, for a given cluster, the best fit red sequence, 
a model red sequence for comparison and the dwarf/giant 
magnitude limits at the corresponding redshift. The model 
red sequence is constructed assuming that the correlation 
between galaxy colour and magnitude is driven solely by an 
increasing metallicity with stellar mass. We use the Chariot 
& Bruzual (2007) synthetic stellar population code (here- 
after CB07) with a Salpeter initial mass function and a single 
burst of star formation at z — 3 to compute synthetic galaxy 
colours over a range of metallicities. Each model CMR also 
contains a small contribution (3% by mass) of ultra metal- 
poor stars, which Maraston et al. (2009) find is necessary to 
reproduce the (g' — r') colour of early type galaxies at similar 
redshifts. Without this component the models are up to 0.1 
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magnitudes too blue. The model is then calibrated to repro- 
duce the magnitude vs metallicity relation Kodama & Ari- 
moto (1997). It can be seen from Figure |A"T1 that the model 
CMR is in qualitative agreement with the observations, in- 
cluding their evolution. The model CMR slope however is 
typically shallower than the best fit CMR and this difference 
in slope can vary appreciably from cluster to cluster. This 
is likely a consequence of the over-simplicity of calibration 
to a fixed magnitude metallicity relation and the use of a 
single stellar population. We take the model as illustrative 
only and for the purposes of classifying giants and dwarfs 
we use the CMR as defined by the fit. 
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Figure Al. Illustration of the red sequence fitting procedure for 20 clusters selected randomly from our sample. They are arranged in 
increasing rcdshift from the top left. The solid red line shows the CMR according to a linear fit to the data while the dotted magenta 
line shows the CMR defined by a stellar population that was created in a burst at z = 3 and aged passively (see text). The dashed red 
lines show the colour offset A(</' — r') = ±0.2 used for selecting red sequence members. The vertical blue lines show the limits to classify 
dwarfs and everything brighter than this is classified as a giant. 
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